arXiv: 1509.0205 lv2 [astro-ph.GA] 11 Feb 2016 


Mon. Not. R. Astron. Soc. 000.^(^(0000) Printed 12 February 2016 (MN style file v2.2) 


Gas Inflow and Metallicity Drops in Star-forming Galaxies 

Daniel Ceverino^’^*, Jorge Sanchez Almeida^’'^, Casiana Munoz Tunon^’'^, 
Avishai DekeP, Bruce G. Elmegreen®, Debra M. Elmegreen^, Joel Primack® 

^ Centro de Astrobiologia {CSIC-INTA), Ctra de Torrejon a Ajalvir, km 4, E-28850 Torrejon de Ardoz, Madrid, Spain 
^Astro-UAM, Universidad Autonoma de Madrid, Unidad Asociada CSIC, E-28049 Madrid, Spain 
^Instituto de Astrofi'sica de Canarias, E-38205 La Laguna, Tenerife, Spain 
^Departamento de Astroflsica, Universidad de La Laguna, E-38205 La Laguna, Tenerife, Spain 

^ Center for Astrophysics and Planetary Science, Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 
^IBM Research Division, T.J. Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, NY 10598, USA 
"^Department of Physics and Astronomy, Vassar College, Poughkeepsie, NY 12604, USA 
^Department of Physics, University of California, Santa Cruz, CA, 95064, USA 


ABSTRACT 

Gas inflow feeds galaxies with low metallicity gas from the cosmic web, sustaining 
star formation across the Hubble time. We make a connection between these inflows and 
metallicity inhomogeneities in star-forming galaxies, by using synthetic narrow-band images 
of the Ha emission line from zoom-in AMR cosmological simulations of galaxies with stellar 
masses of M* ~ 10® Mq at redshifts z = 2 — 7. In ^50% of the cases at redshifts lower 
than 4, the gas inflow gives rise to star-forming, Ha-bright, off-centre clumps. Most of these 
clumps have gas metallicities, weighted by Ha luminosity, lower than the metallicity in 
the surrounding interstellar medium by ^0.3 dex, consistent with observations of chemical 
inhomogeneities at high and low redshifts. Due to metal mixing by shear and turbulence, 
these metallicity drops are dissolved in a few disc dynamical times. Therefore, they can be 
considered as evidence for rapid gas accretion coming from cosmological inflow of pristine 
gas. 
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1 INTRODUCTION 


A number of observational works have detected chemical 
inhomogeneities i n the gas of star-f o rming galaxie s , both 
at high redshift (ICresci et al.l 120101: lOuevrd et al.l 12012h 
and in the local universe (Richards et al. 2014; see 
Sanchez Almeida et al. 2014a for a review). Most 
galaxies develop negative metallic ity gradients, with the 
metal licity decreasing inside-out (Edmunds & Greenhowl 

1995t __ Ivan Zee et al. _ _ 1998 : 

Magrini et akl 120071: ISanchez et al.l 120141 IStott et aP 2014; 


Kehrig et al.l l2015h . On top of such large-scale variation. 


many galaxies show localized, kpc-size starbursts with 
the metallicity in the burst lower than in the surrounding 
interstellar medium (ISanchez Almeida et al.l 120131 1201 4bl 
1201 5ll . These metallicity drops are of the order of 0.3 dex or 
larger. Since the time-scale for gas mixing in disk galaxies 
is expecte d to be short, on the order of the rotational period 
or shorter ([^ Avillez & Mac Lo^l2002t lYang & Krumholzl 
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1201 21: iPetit et akl 1201 5l) . the low-metallicity gas sustaining 
the starbursts must have arrived to the galaxies recently. 

A plausible interpretation of these observations is the 
triggering of sta r-formation by gas recently accreted from 
the co s mic w eb (ICresci et al.ll2010l : ISanchez Almeida et al.1 
1201 31 1201 5ll . Numerical simulations predict that the 
inflow of cosmic-web gas has sust ained star formati on in 
galaxies along the Hubble time (iDekel et alJ 1201 3l and 
references therein). The gas being accreted is metal-poor 
(Ivan de Voort & Schavel 12012) and it is e xpected to induce 
giant star-forming clumps in early disks dPekel et al.ll2009l : 
ICeverino et al.ll201^ iMandelkeret al.ll2014t) . The ms from 
the web may trigger a starburst as it gets compressed when 
approaching the disk, or it builds up the gas mass in the 
disk which eventually gives rise to starbursts due to internal 
instabilities. Only if the buildup and triggering occurs rather 
fast, in a time-scale similar to or shorter than the mixing 
time-scale, the recently accreted, low-metallicity gas can 
still dominate the mass of the star-forming clumps. 

This paper studies whether this conjecture is valid 
in cosmological simulations of galaxies, namely whether 
metal-poor star-forming gas clumps in high-z protogalaxies 
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can be associated with individual gas accretion events. 
IVerbeke et al.l(l2014t) simulate the effect of gas clouds falling 
into disks of various masses and types. Their conclusion 
is that gas accretion may give rise to low metallicity 
starbursts depending on details of the accretion event, such 
as the accreted gas mass, the impact parameter of the 
stream, or whether the orbit is prograde or retrograde. These 
simulations are idealized with ad-hoc initial conditions. 
We want to see whether the phenomena of metallicity 
drops associated with cosmic gas accretion appear in a 
cosmological setup, where gas infall is not imposed but 
appears naturally, as a result of the formation of large-scale 
structures. 

We aim at studying and characterizing the physical 
process, rather than at explaining a particular set of 
observations. However, we want to make predictions that 
connect a given observable such as the gas metallicity 
weighted by Ha luminosity, with cosmological gas inflow. 
In section ij2] we describe the sample of simulations. Sj3] 
describes the mock Ha images of the sample. |j4] discusses 
the metallicity drops found in the simulations and we 
connect them with recent gas accretion events in 0 Finally, 
0is devoted to discussion and summary. 


Table 1. Properties of the zoom-in simulations, ordered by decreasing 
gas fraction. Columns show name of the run, redshift, virial radius, virial 
mass, stellar mass, gas mass within 0.i/?vir SFR in Mq yr“^. 
Masses are in Mq and radii in kpc. Virial propert ies are computed 
using an overdensity, A (z)2:^18Q, above mean density jPekel et alj|2013l : 
iBrvan & NormaiJl998h . 


Run 

Z 

R-v'ir 

virial 

mass 

Stellar mass 

gas mass 

SFR 

V05 

3.0 

31 

5.3 

X 

lOio 

4.0 

X 

10® 

5.9 

X 

10® 

0.5 

V09 

5.2 

20 

5.2 

X 

lOio 

5.2 

X 

10® 

6.0 

X 

10® 

2.2 

V22 

6.1 

17 

5.1 

X 

lOio 

6.2 

X 

10® 

6.8 

X 

10® 

4.3 

V07 

4.9 

18 

3.3 

X 

lOio 

5.2 

X 

10® 

5.7 

X 

10® 

3.0 

VOl 

3.5 

26 

4.5 

X 

lOio 

4.5 

X 

10® 

4.6 

X 

10® 

1.2 

V04 

1.9 

55 

12 

X 

lOio 

10 

X 

10® 

10 

X 

10® 

2.6 

V06 

5.2 

20 

5.5 

X 

lOio 

5.9 

X 

10® 

5.8 

X 

10® 

3.2 

V26 

4.9 

22 

5.9 

X 

lOio 

5.7 

X 

10® 

5.3 

X 

10® 

4.0 

V17 

6.7 

16 

5.1 

X 

lOio 

5.6 

X 

10® 

5.0 

X 

10® 

5.0 

V13 

3.3 

33 

8.1 

X 

lOio 

7.8 

X 

10® 

6.9 

X 

10® 

2.3 

V02 

2.7 

37 

7.0 

X 

lOio 

8.1 

X 

10® 

5.9 

X 

10® 

1.1 

V25 

3.5 

32 

8.8 

X 

lOio 

8.3 

X 

10® 

5.7 

X 

10® 

2.2 

V23 

4.0 

27 

6.6 

X 

lOio 

10 

X 

10® 

6.2 

X 

10® 

2.8 

V27 

3.8 

32 

9.9 

X 

lOio 

8.9 

X 

10® 

5.1 

X 

10® 

0.9 

VOS 

3.2 

33 

7.1 

X 

lOio 

9.9 

X 

10® 

4.8 

X 

10® 

1.4 

V33 

4.5 

25 

7.6 

X 

lOio 

9.5 

X 

10® 

4.5 

X 

10® 

2.4 

V20 

6.1 

16 

3.9 

X 

lOio 

8.3 

X 

10® 

3.5 

X 

10® 

2.3 


2 SIMULATIONS 


The zoom-in simulations used in this p aper were drawn from 
a larger dataset, th e Vela sample (ICeverino et al.l 1^14 
IZolotov et al.ll20L^ . We select for each simulation the latest 
snapshot that fulfills the two following criteria: 

• A galaxy stellar mass of M* < 10® Mq, 

• A specific star-formation-rate of sSFR ^ 1 Gyr“^, 


where sSFR = SFR/M^, (Table [T]i- Therefore, we select 
star-forming galaxies with low masses (M, = (0.4 — 
1) 10® Mq), and gas rich (gas-over-star mass ratio about 
unity) at high-redshifts {z = 6.7 — 1.9). These protogalaxies 
are expected to undergo large, quasi-continuous gas inflows, 
so that they are the best candidates for setting the connection 
between cosmological gas inflow and metallicity drops. 
Ongoing or recent major mergers with strong inflows of stars 
and dark matter could complicate this connection, so that 
they are excluded from the analysis. 

The simulations were performed w ith the ART code 
jKravtsov et al.l Il997t I Kravtsov! l2003h . which accurately 
follows the evolution of a gravitating N-body system and 
the Eulerian gas dynamics using an AMR approach. Beyond 
gravity and hydrodynamics, the code incorporates many 
of the physical processes relevant for galaxy formation. 
These processes, representing subgrid physics, include 
gas cooling by atomic hydrogen and helium, metal and 
molecular hydrogen cooling, photoionization heating by 
a constant cosmological UV background with partial 
sel f-shielding, star formation and feedback, a s desc ribed 


in 


.Ceverino & Klvpin (l2009h . ICeverino et al.l (l2010t) . and 
ICeverino et al.l (12014 ). In addition to thermal-energy 
feedback, the simulations use radiative feedback. This model 
adds a non-thermal pressure, radiation pressure, to the 
total gas pressure in regions where ionizing photons from 
massive stars are produced and t rapped. In the curren t 
implementation, named RadPre in ICeverino et alJ (l2014t) . 


radiation pressure is included in the cells (and their closest 
neighbors) that contain stellar particles younger than 5 Myr 
and whose gas column density exceeds 10®^ cm“®. 

Metals are pr oduced and adv ected according to the 
model described in iKravtso^ (l2003l) . The code follows two 
kinds of metals: metals produced in SN-II (Oxygen and 
other a-elements), and metals produced in SN-Ia (Iron-like 
elements). Stars with masses between m* = 8 Mq and 
100 Mq eject /z m* mass of metals in core-collapsed 
supernovae, where /z = min(0.2,0.01 — 0.06), which 

crudel y approximates the SN yields of IWooslev & Weaverl 
(1 19951) . After integrating a Chabrier IMF, this translates into 
the ejection of ^ 1% of the initial stellar mass of a particle 
representing a single stellar population. These metals are 
injected in the cell hosting the star particle at a constant rate 
during 40 Myr, the time that takes the last SN-II to explode. 

The initial conditions of these runs contain from 6.4 
to 46 10® dark matter particles with a minimum mass 
of 8.3 X 10"^ Mq, while the particles representing single 
stellar populations that were formed in the simulation have a 
minimum mass of 10® Mq. The maximum spatial resolution 
is between 17-35 pro per pc. More deta i ls can be found in 
ICeverino et al.l (l2014l) and IZolotov et al.l (l2015l) . 


3 Ha IMAGES 

If we want to make a connection between the gas inflow and 
the metallicity drops in starburst regions, the first step is to 
generate synthetic narrowband observations of Ha emission 
that highlights star-forming regions. The mock Ha images 
of the sample are genera ted according to the method outlined 
in ICeverino et al.l (l2012l) . In short, we first compute the star 
formation rate density, psFR, using the distribution of stellar 
particles younger than 50 Myr. Then, we compute the Ha 
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V09 V27 

Figure 1. Synthetic Ho maps of galaxies with off-centre bright clumps, marked by a white circle with a radius of 0.4 kpc. A white dot mai'ks the galaxy center 
of mass. The size of each image is 10 x 10 kpc^. In most of these runs, the clumps have a metallicity lower than in the inter-clump medium (red labels). A 
single case with a slightly off-center clump with metallicity higher than in the inter-clump is labeled in black. 



Figure 2. Synthetic Ho; maps of galaxies with only one bright central clump marked by the white circle. The size of the images and the labels are as in 
Figure^ 2 out of 11 cases show metallicity drops at the galaxy centre. 
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emissivity e Hr>, based on the iKennicuti (Il998h conversion, 
adjusted to a lChabrie^ (l2003h IMF, 

logCHa = log PSFR +41.33, (1) 

where CRa is in ergs“^kpc“^ and psFR is 
in Mq yr“^ kpc“^. This equation holds as long as the Ha 
photons trace the underlying star-formation events, which is 
the ca se on scales of a few hundred parsecs dKennicutt et al.l 
l2007h . The Ha surface brightness is obtained by integrating 
the emissivity along a line-of-sight parallel to the galaxy 
rotation axis (face-on view). Radiative transfer effects and 
dust extinction are not considered. However, these effects 
are small in th ese dwarf galaxies , where the dust column 
density is low (iFisher et al.ll2014l Sanchez Almeida et al., 
submitted). A Gaussian smoothing of FWHM=0.3 kpc is 
implemented in the hnal images. 

The galaxies are classihed in two groups, according 
to their morphology in the face-on, Ha images. The hrst 
group looks clumpy in Ha, showing one or more Ha-bright, 
off-centre clumps (Figure[T]). The second group shows only 
one Ha-bright clump at the galaxy centre (Figure|2l). These 
Ha-bright clumps could be associated with the starburst 
regions observ ed in star-forming galaxies at high and 
low re dshifts (iGenzel et alJ 1201 It ISanchez Almeida et akl 
l2014iJ) . 

Around one-third of the sample of star-forming galaxies 
are clumpy. This clumpy fraction increases to ^50% if we 
exclude redshifts higher than 4. This is consisten t with UV 
obser vations of bright clumps in high-z galaxies (iGuo et akl 
I 2 OI 5 I) . Their hgure 10 shows a 60% fraction of clumpy 
galaxies with stellar masses of log(M*/M0)=9-lO in the 
redshift range z = 0.5 — 3, slightly higher mass es than in our 
sample , (see also lElmegreen et al.l (l2007h and iTadaki et al.l 

(Eoil). 


4 METALLICITY DROPS 

The next step is to hnd the metallicity drops by computing 
the mean gas metallicity, weighted by Ha emissivity, of the 
Ha-bright clumps and the inter-clump medium. We translate 
the gas density of metals produced by supernova type-II 
into oxygen abundance s by using the formula described in 
iMandelker et akl (l2014l) . 

^ _ /o ZSNII _ 0-5 Z^SNII - 2 ') 

'H~ XA “ 0.755 x 16 ^ 

where zsnii is the metal-to-gas mass ratio within a cell, 
X = 0.755 is the hydrogen mass fraction, and A = 16 
is the atomic weight of oxygen. We also assume that half 
of the mass in metals released in SNII are oxygen atom s 
(fo = 0.5), roughly based on IWooslev & WeavCT (Il995h . 
The actual values of the scaling factors do not affect our 
results which are based on metallicity ratios rather than on 
the absolute value of the metallicities. 

The Ha-weighted gas metallicities are measured in 
circular apertures of radius 0.4 kpc, centered on the brightest 
off-center clump in Ha for each clumpy galaxy (Figure [T]i. 
For the non-clumpy galaxies (Figure [2]i, the aperture is 
centered at the Ha-bright clump at the galaxy center of 
mass. In this first approach, we ignore all the smaller 
(dimmer) clumps, because they are usually not observed 
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Figure 3. Metallicity versus SFR surface density measured in apertures 
of 0.4 kpc radius. Each galaxy has two measures, linked by a dotted line, 
labeled by the iTin ID. The highest SFR value in each galaxy corresponds to 
the Ha-bright clump. The lowest value corresponds to a randomly placed 
aperture. The small dots represent 100 random apertures for V02 (top) and 
V13 (bottom), excluding the Ha-bright clump and the galaxy centre. The 
squares are joined by red or black lines depending on whether the clump 
is metal poorer or metal richer than the inter-clump medium. Cyan squares 
and lines correspond to the center of clumpy galaxies. 
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in the galaxies of this mass scale. Therefore, we focus on 
the brightest clump that us ually corresponds to the head in 
observed tadpole gala xies (lElmegreen & Elmegreenll2010t 
lElmegreen et al.ll2012h . This is where a large amount of gas 
is needed to feed high SERs. 
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Figure 4. Face-on maps of gas surface density (E), velocity in the X-axis of the galaxy plane (Vx), metallicity (O/H), velocity in the Y-axis (Yy), gas 
metallicity weighted by Ho emissivity (O/ HjfP, and radial inflow of V27 at ^ = 3.8. In each panel, three squai'e apertures mark the galaxy center and two 
regions with low metallicity at the end of two incoming streams of pristine gas. Both apertures are close to Ho-bright star-forming regions. The size of the 
images is 20 X 20 kpc^. 


In addition to the aperture centered on the clump, 
we also consider another aperture randomly placed in the 
inter-clump medium with a SFR surface density close to 
the galaxy median value. Therefore, each galaxy has two 
metallicity measures. The apertures centered on Ha-bright 
clumps have a high SFR surface density, much higher 
than in the inter-clump medium. The actual position of 
the aperture chosen for reference does not change the 
conclusions significantly, as shown by 100 different random 
apertures for V02 and VI3. 

The top panel of Figure |3 shows the metallicity in the 
clumpy galaxies of Figure [T] The galaxies with metallicity 
drops are shown using red symbols and lines in Figure [3] 
In almost all cases, the gas metallicity in the clump is 
significantly lower than the metallicity of the inter-clump 
medium. The median drop in metallicity is 0.33 dex with 
a maximum value of 0.5 dex (V02). The only case of an 
off-centered clump without a metallicity drop is shown in 
VOS. This clump is also the nearest to the center of mass, so 
this could be an intermediate case between off-centered and 
central clumps. 

The central clumps of Figure |2] have higher SFR 
densities than the off-centered clumps, most probably 


because these galaxies are more compact than their clumpy 
counterparts. The metallicity of the central clump is similar 
to or higher than the metallicity of the galaxy. More 
than 50% of the random apertures described above have 
metallicities lower than the central value. This is consistent 
with the fact that irregular galaxies with stee p inner profiles 
could have negative metallicity gradients (iPilvugin et al.l 
IMsh . VI and V7 are two exceptional cases that show 
a central metallicity which is ^ 0.5 dex lower than the 
metallicity at larger radii. In these cases, large radial flows 
of gas decrease the metallicity at the center. 


It is remarkable that the Ha-bright, off-centre clumps 
have lower metallicities than other star-forming regions with 
lower Ha luminosities, because these clumps have high 
star formation rates and they are regions where significant 
amount of metals are injected into the gas due to supernova 
explosions. These regions must have accreted large amounts 
of gas with low-metallicity, possibly coming from streams 
of inflowing gas. This inflow of low-metallicity gas could 
explain these inhomogeneities in metallicity. 
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Figure 5. Evolution of the gas surface density of V27, in time intervals of 10 Myr. Time increases from top-left to bottom-right. The rightmost panel of the 
third row (mai'ked with the direction of the streams) coincides with the top-left panel of Figure |4] The three black squares are located at the same positions as 
in Figure |4] Dense, low-metallicity stream gas is coming from the right side of the view (marked by an ellipse) and accumulates material in the top-left side 
of the galaxy, increasing the gas surface density and the SFR in that region. These star-forming clumps are dispersed after '--'50 Myr, i.e., a few times the disc 
dynamical time. 


5 GAS INFLOW 

The inflowing gaseous streams join the disk through an outer 
ring, that can be tilted at large radii and then gradually 
aligne d with the galaxy disk at low radii (iDanovich et akl 
l2015h . This implies that the streams join the disk generally 
within a plane and at a given point in angle (j), which may 
change as the streams-plus-galaxy system evolves. Figure |4] 
shows an example of gas accretion through streams (V27 


at z = 3.8). Two streams coming from the circumgalactic 
medium can be seen at the left and right sides of the maps. 
They are characterized by their low metallicities, high speed 
and inflow. The streams meet the galaxy around the regions 
marked with two off-centered squares in Figure |4] In these 
regions, gas accumulates and drives high SFR and high 
Ha luminosities (Figure [I}. The Ha-weighted metallicity 
of these Ha-bright clumps is much lower than the average 
metallicity in the galaxy. Therefore, these metallicity drops 
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Figure 6. Face-on maps of gas surface density (S), mass-weighted gas 
metallicity {OJH), gas metallicity weighted by Ha emissivity {O jHuP, 
and radial inflow of V02 at 2 = 2.7. In each panel, two square apertures 
mark the galaxy center and the off-centered clump. The size of the images 
is 20 X 20 kpc^. The color palette is the same as in Figure|4] 


mark regions where low-metallicity streams accumulate 
material in the galaxy. 

Figure |5] shows the temporal evolution of the gas 
in time steps of 10 Myr, 110 Myr before the snapshot 
shown in Figure H) In this sequence, relatively dense 
low-metallicity gas is coming from the right side of the 
view. This stream of inflowing gas accumulates material in 
the top-left side of the galaxy, increasing the gas surface 
density and the SFR in that region. Therefore, gas accretion 
into galaxies promotes an Hcr-bright burst of low-metallicity. 
After ~50 Myr, a few times the disc dynamical time, 
these star-forming clum ps dissolve due to t he combination 
of radiati ve feedback (IMoodv et al.l I20l4l and turbulent 
motions (lYang & Krumholzl l2012h . The low metallicity 
drops are signatures or hints of recent gas accretion into 
galaxies. 

Figure |6] shows another example of the link between 
metallicity drops and gas inflow. In V02, the off-centered 
clump is much closer to the galaxy centre than in the 
previous case. As shown in Figure @1 low-metallicity gas 
wraps around the galaxy. The metal-poor clump, according 
to its Ha-weighted metallicity, is located at the end of one 
of the inflowing streams (at the bottom of the clump). Other 
cases show similar patterns of gas accretion. 

We can quantify this gas accretion by following the 
evolution of the gas mass within a sphere of fixed radius, 
0.15i?vir- Figure I 2 ] shows this temporal evolution for gas 
and stars of the clumpy and non-clumpy galaxies shown in 
Figure [T] and Figure |2l Masses are normalized to the value 
at the selected snapshot (Table[T]i. We choose a time interval 
of ±0.1 Gyr around the snapshot, so we can quantify the 
inflow (or outflow) before (and after) that moment captured 



-0.1 -0.05 0 0.05 0.1 

t-to (Gyr) 


Figure 7. Evolution of the gas mass (top) and stellar mass (bottom) 
within 0.15/?vir of the galaxies shown in Figure ^ and Figure [2] In all 
the five clumpy galaxies with metallicity drops induced by gas accretion 
(labeled), there is a significant increase in the amount of gas within ~0.1 
Gyr before the selected snapshots at £o- Cyan dotted lines correspond to 
the non-clumpy galaxies. There is no significant stellar accretion with the 
exception of the mergers between 0.05 and 0.09 Gyr after £o- 

in Figure [T] and we allow some time for the formation of 
clumps from this inflow material. 

In three out of five cases with metallicity drops, there is 
a steep increase in the gas mass in a period of 0.1 Gyr before 
the selected snapshot. The gas mass doubles within this 
period. The increase in gas is not related with merger events, 
because there is no sudden increase of the stellar mass within 
the same period. Instead, the stellar mass increases slowly. 
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about 10-20%, due to star-formation within the selected 
sphere. A exception is a merger event in V27 at 0.08 Gyr. 
Therefore, the increase of gas mass in these galaxies with 
metallicity drops is due to the continuous infall of gas from 
cosmic streams. 


6 SUMMARY AND DISCUSSION 


We addressed the continuous gas accretion of low metallicity 
gas into star-forming galaxies at different redshifts, using 
a set of zoom-in AMR cosmological simulations of 
galaxy formation. In ^50% of the cases at redshifts 
lower than 4, the accretion gives rise to star-forming, 
Ha-bright, off-centre clumps. Most of these clumps have 
gas metallicities, weighted by Ha luminosity, lower than the 
average metallicity in the surrounding interstellar medium. 
The typical values of these metallicity drops are around 
0.3-0.5 dex, consistent with observations d&esci et al.ll2010l: 
ISanchez Almeida et al.l 12013112014bi l2015h . Therefore, the 
observed metallicity drops can be considered as evidence 
for rapid gas accretion coming from cosmological inflow of 
low-metallicity gas. 

Low metallicity inhomogeneities are dispersed in 
^^50 Myr, only a few times the disc dynamical time, 
td = Td/vd = 2 kpc/100 kms“^ = 20 Myr, where 
rd and Vd are the characteristic radius and rotational 
velocity of these galaxies. Therefore, low-metallicity gas 
gets mixed with galactic interstellar medium in less than 
an orbital period. These short time-scales are consistent 
with the time-scales for dispersio n of non-axisym metric 
inhomogeneities in galaxy discs (iPetit et al.l l2015h . The 
metal mixing is led by shear and turbulence, driven 
by gravitational disc instabilities. Shear accelerates the 
diffusion of metals by transferring inhomogeneities from 
large scales to small scales, where tur bulence is able to 
dissip ate them in a dynamical time-scale (lYang & Krumhol3 
Hall). Feedback may also contribute to the dispersion 
of met als by the disrup tion of the small metal-poor 
clumps (iMoodv et al.ll2014l). More massive galaxies hosting 
giant, long-lived clumps (iDekel e t al.ll200 9t ICeverino et al.l 


l2010t iMmdelker et alJ l2014t iBournaud et al. 2014t) . may 
have metallicity drops with longer lifetimes, if clump 
self-enrichment due to supernovae yields is not a relevant 
process. Future work focused on these giant clumps will 
clarify these issues. 
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